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We study some effects arising from periodic modulation of the asymmetry 

and the barrier height of a two-well potential containing a Bose-Einstein con- 

^D ' densate. At certain modulation frequencies the system exhibits resonances, 

" 

^^ ' which may lead to enhancement of the tunneling rate between the wells and 

1 ^ ' which can be used to control the particle distribution among the wells. Some 

of the effects predicted for a two- well system can be carried over to the case 

of a Bose-Einstein condensate in an optical lattice. 
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^ ! Since the experimental realization of Bose-Einstein condensates (BEG) one has consid 

Oh! ered the possibility of extending one-mode models to two or more modes [1]. This raises the 
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I. INTRODUCTION 



issue of the relative phase between modes. As indicated by several authors [2], a two- well 



Q^, BEG may exhibit features that are not covered by the semiclassical description in terms of 



the Gross-Pitaevski equation. These features are significant at low particle numbers and 
for strong interactions. In previous work, we discussed some aspects of the dynamics of a 
two- well BEG in the strong-interaction regime [3] . This is close to the experimental situation 
for a BEG in a double- well trap, designed in Ref. [4]. 

A sensitive way to probe the properties of a BEG in a double-well potential with strong 
interatomic interactions is to look for resonant behavior when a parameter of the system is 
periodically modulated. The response of the system may be expected to be very sensitive 
to the value of the modulation frequency in the neighborhood of a resonance. A periodic 
perturbation can be implemented in various ways. One example would be periodically 
modulating the trapping potential. Salmond et al [5] study a numerical model of a double- 
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well potential with periodically modulated coupling between the wells. This semiclassical 
analysis reveals the existence of uncoupled regions with chaotic and regular motion. The 
inclusion of the quantum nature of the evolution leads to transitions between these regions. 
Another type of periodic perturbation can be imposed by periodically modulating the energy 
difference between the ground states in the two wells. 

Periodic modulations are known to give rise to dynamical localization in some cases. This 
effect has been widely discussed in the literature in the case a particle in a periodic potential, 
such as an electron in a crystal or an atom in an optical lattice [6] and [7]. When the particle 
also feels a uniform force in addition to the lattice potential (a tilted optical lattice), it is 
found to have an infinite discrete set of equidistant energy levels, with a level separation 
that is determined by the strength of the uniform force [7] . A variation of the magnitude 
of the uniform force affects the phase of the state. So, when this magnitude is periodically 
modulated, resonances may be expected. The population in one well is described by adding 
the amplitudes for arriving at this well from various other wells, each one with a different 
phase shift. In the long time limit, when the time of observation is long compared to the 
period of the modulation, this gives rise to destructive interference, leading to a suppression 
of the net tunneling rate. Hence, the asymptotic distribution over the wells coincides with 
the initial one, and dynamical localization has been realized. 

Another example of dynamical localization arises for a single two-level atom in a driv- 
ing field with a periodically modulated detuning [8]. For certain ratios of the modulation 
frequency and the strength of the field, the atom is localized in its ground state. The time 
scale is restricted by the modulation frequency. 

In the present case of a BEC in a two-well potential with a fixed total number of atoms, 
the state space is finite dimensional. In the SU(2) representation of the operator algebra, 
the Hamiltonian has a quadratic term due to the interatomic interactions. So, in this sense 
the system is quite different from that of an atom in a tilted lattice, with its infinite number 
of states and a Hamiltonian that is linear in the SU(2) operators. Still, there are some 
obvious similarities: the discrete structure of the energy and the presence of interatomic 



interactions and tunneling between wells as competing processes. Therefore, we expect 
interesting effects also in the two- well case when the energy difference or the hopping between 
wells is periodically modulated. From a formal point of view, the analysis in the present 
paper may be regarded as a generalization of the process of dynamical localization for the 
Hamiltonian with a quadratic term. Specifically, this paper considers the possibilites of 
coherent control of a BEC in a double-well potential by using any kind of time-periodic 
perturbation. 

II. BEC IN A DOUBLE POTENTIAL WELL 

We describe a BEC in a double potential well in terms of a one-particle Hamiltonian H^^^ 
and a two-particle interaction U{f,f'). The states |1) and |2) are the the localized ground 
states in either well, with wave functions ipi{r) and ip2{'f^- On the basis of the states |1) and 
1 2), the one-particle Hamiltonian has the matrix elements 

(1| ijW |1) = - (2| i/W |2) = he/2, (1| ifW |2) = (2| i/^ |1) = -M/2 . (1) 

In the case that e = 0, the coupling between the wells lifts their degeneracy, and creates an 
energy sphtting h6 between the even ground state \g) and the odd excited state |e), defined 
by 

l^) = ^(|l) + |2)); |e) = i=(|l)-|2)). (2) 

When we restrict ourselves to these two states, the field operator in second quantization has 
the standard form 

^{T')=aiiJi{r)+a2Mr), (3) 

with tti the one-particle annihilation operator in the two states, which together with the 
corresponding creation operators obey the bosonic commutation rules. When we substitute 
this expression in the formal expression 



H = f rfr$1"(r)if ^§(7^) + - I d-rd-r'^\-r')^\Y")U{f, r')§(r)§(r'). (4) 

for the second-quantized Hamiltonian, we find 

H = '^h{i\ H^'^'^ \k) a}{dk + - ^ h{i,k\U \l, m) Ki^k,i,malalaiam, (5) 

i,k i,k,l,m 

where the indices i, j, k, I = 1 or 2, and the matrix elements are taken between the states 
V'l and i/j2- ■ 

At sufficiently low energy, the two particle interaction is well approximated by the contact 
potential U{f,r') = {A7Th'^a/m)6{f — f), with a the scattering length. The function ipi and 
■?/'2 states have the same form, and we assume that they do not overlap. So we obtain the 
following expression for the Hamiltonian 

-- he /^+^ ^1^ \ h6 /^f^ ^t^ \ hK, /^f-t-- -^ ^i^i^ ^ \ ,„. 

H = — I a[ai — 0202 I I a,ia2 + 02^1 1 H I a[a|aiai + a2'22'^2fl2 1 • (oj 

where the parameter k, defined by 

h^ = ^^ [df\Mr)\^ (7) 

171 J 

measures the strength of the interatomic interaction. 

For convenience we express the Hamiltonian (6) in terms of SU(2) operators by applying 
the standard Schwinger representation of two modes. This leads to the definition 

Jz = - ( ajoi — 02^2 ] 1 J+ = a|a2 , J~ = alai. (8) 

These operators are related to the Cartesian components of a fictitious angular momentum 
by the standard relations J± = Jx i: iJy They obey the commutation rules for angular 
momentum operators 

[Z, J±] = ± J± , [J+, J_] = 2Z , (9) 

which generate the SU(2) algebra. These operators commute with the operator for the total 
number of particles A^ = ajai + 0302. The Hamiltonian (6) can be rewritten in the form 



H = H^ + ^(^N'-2n) , (10) 

where the A^-particle Hamiltonian Hn is defined by 

Hn = heJz - hSJ^ + HkJI. (11) 

For a given number of particles A^, the last term in Eq. (10) is a constant, and it suffices 
to consider the dynamics of the subspace of the A^ + 1 number states |n, A^ — n), with 
n = 0, 1, . . . A^, with n particles in well 1, and N — n particles in well 2. This subspace has 
the structure of the angular momentum states, with J = N/2, and the 2J + 1 magnetic 
quantum numbers fi = n — N/2, with /i = — J, — J+1,...,J. Note that /i is half the 
difference of the particle number in two wells. For a given particle number A^ we represent 
the number states by the quantum number /x, so that \fi) = \n, N — n) . The action of the 
operators Jq and J± on the Fock states has the well-known behavior 

j,|/i) = /i|/i) , j±|/i) = v^(j?70G7±77TTy|/i±i). (12) 

This also determines the action of the Cartesian operators Jx and Jy. 

III. QUANTUM STATES IN TWO WELLS 

The Schwinger representation of the operators occurring in the Hamiltonian suggests in 
a natural way various possible choices of states of A^ atoms in the two wells. Arecchi et al [9] 
introduced the spin coherent states (SCS) [10], in analogy to the Glauber coherent state of 
a mode of the quantum radiation field. The SCS follow from applying an arbitrary rotation 
to the state |/i) with n = J. As rotation operator we take 

R{9,(f)) = exp{—i(l)Jz) exp{—i9Jy) ex\i{i(f)Jz) = exp[— i9{Jy cos (j) — J^^sin^)], (13) 

which represents a rotation over an angle 6, around an axis in the xy-plane, specified by the 
angle with the y-axis. The SCS \9, 0; J) is 

1^,0; J) = %, 0)1 J) (14) 



which is also the eigenstate with eigenvalue J of the component "z? ■ J of the angular- 
momentum vector in the direction it specified by the polar angle 9 and the azimuthal 
angle 0. Just as the Glauber coherent states of a mode with annihilation operator c can 
be obtained by acting with a displacement operator D [Q = exp ((^c^ — (*'cj on the vacuum 
state, the SCS follows by a rotation R{9,(f)) = exp (C-^- ~ C*'^+) with ( = (^/2) exp (i0), 
acting on the state \J). When we view this state \J) as the ground state, the operator J+ 
is analogous to the annihilation operator, since J+\J) = 0. An essential difference between 
the two cases is, of course, that the state space of a radiation mode has infinite dimensions, 
while the dimension of the angular- momentum state space is 2J + 1. 

In our case, the analogy is carried one step further, since the SCS defined by (14) do not 
represent angular-momentum states, but refer to the states of A^ atoms, distributed over two 
potential wells. The ground state | J) represents the state with all particles in the first well. 
When we substitute the identity \J) = ( a{ ) \vac) /y N\ with N = 2J into the right-hand 
side of (14), we obtain an expression for the SCS in the language of the two wells, in the 
form 

1/9 9 \^ 

\9, (b\ J) = —== I cos - a] + e^"^ sin - al | \vac) . (15) 

VM V 2 ^ 27'^ ^ ' 

We can interpret (15) as a state with A^ atoms in the one-particle superposition state 

9 9 

cos- 11) + 6^"^ sin- 12) 

of the two wells. A special case arises for 9 = 7r/2, when the average populations of the 
two wells are the same. Then the state (15) describes a collective mode from two interfering 
sources of equal intensity, and its expansion in number states is given by 



I^A0;'^) = ^E 



N I nr \ 



2N/2 

n=0 



N 



AN-n^\n^N-n) (16) 



Such a state can be considered as a state with a well-defined phase difference 0. The atom 
distribution over the two wells is binomial, and they have been termed phase states (PS) of 
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a two- mode boson system in Ref. [11]. For simplicity, we suppress the value 7r/2 in this case, 
and we simply denote the PS as |0; J) . Upon rotation around the z-axis, a PS transforms 

as 

exp{—iaJz) 10; J) = exp{—iaJ) |0 + a; J) (17) 

The concept of Glauber coherent states of a radiation mode has been generalized by 
de Oliveira et al [12], who introduced so called displaced coherent states defining them 
as a displaced number state, rather than a displaced vacuum state. The corresponding 
generalization of a SCS is found when the rotation operator (13) acts on a number state \fi). 
The resulting displaced spin coherent states (DSCS) are 

\9,<P;fi)=R{9,<p)\fi). (18) 

They are the eigenstates of the angular-momentum component it ■ J with eigenvalue fx. In 
the special case that 6 = 7i/2 and = 0, we find that R{6,(j)) = exp(— ^7^/2), and this 
component is simply J^. Its eigenstates are denoted as 

and they obey the eigenvalue relation J^ |/i)^ = /i l/i)^,. The state |/i)^ describes a state with 
J + fi atoms in the even state \g), and J — fi atoms in the odd state |e). These states are 
coupled by the ladder operators 

J^ = exp(— ■J7rJy/2)J±exp(i7rJy/2) = — J^ ± iJy, (19) 

according to the relations 

When 9 = n/2 and = 7r/2, the DSCS are indicated as 

|7r/2,7r/2;/i) = |/i)^. 



which are eigenstates oi it ■ J = Jy, as specified by the relation Jy |/i) = n |/i) . The 
corresponding ladder operators are 

Jy = exp(i7r J^/2) J± exp{-inJ^/2) = 4 =F iJz- (20) 

In the special case that fi = J, the DSCS |/i)^ and |/i) become the SCS | J)^ and | J) , which 
are also the PS with = and = 7r/2, respectively. 

IV. EVOLUTION IN LIMITING CASES 

For a given number A^ of atoms, the evolution is characterized by an evolution operator 
that is governed by the Hamiltonian (11), and that obeys the Schrodinger equation 

ih^ = HnU. (21) 

In order to get an intuitive insight into the evolution, we first consider two extreme cases, 
which are simple to understand. We assume that the two wells have equal energy, so that 
e = 0. If the interatomic interactions are negligible, the quadratic term in (11) can be 
skipped. For a possibly time-dependent coupling strength 6, the evolution operator is 

t/(t) = exp (iT] (t) jL) , (22) 

with f] (t) = J 6{t )dt the area of the coupling pulse. In the language of angular momentum, 
U represent a rotation over an angle —t] around the x-axis. In this case, the states |/i)^ 
are eigenstates of the evolution operator, so that these states acquire only a phase factor 
exp {irj (t) fi) . An initial state in the form of a single number state \fi) state gets rotated 
by the operator (22) and evolves into a superpostion of number states. At the instant that 
1] (t) = 7r/2 an initial number state has evolved into an eigenstate of the operator Jy. 

Conversely, when the interatomic interactions are strong enough on the scale of tunneling, 
the hopping between the wells can get suppressed [3]. Now, a single number state |/x) only 
acquires a phase factor exp(— i/tyU^t). The evolution operator takes the form 



U {t) = 6-'"^^^^ (23) 

which cannot be conceived as a rotation in the angular-momentum space. Since the eigen- 
values of J^ are discrete, the evolution (23) has revivals. First we consider the situation 
that the number of particles A^ is even, so that the eigenvalues /i of J^ are integer. Then 
the eigenvalues of U are exp(— m/i^t) = 1 when t = mT, with m an integer, and T = 2-11 /k, . 
At these times the initial state is reproduced, which proves that the evolution of any initial 
state is time periodic, with period T. For a time t = T/2, which is half the period, the 
eigenvalues of f/(T/2) are exp(— ivryU^) = (— l)'^, which proves that the evolution operator at 
this instant is equal to 

U{T/2) = exp{-iTrZ)- 
For an initial PS |^(0)) = |0; J), we find that the state at the time t = T/2 is 

|^(r/2)) = exp(-mJ)|0;J), 

which is just the opposite PS. At other instants of time, that are a simple rational fraction 
of T, an initial PS can be transformed into a linear combination of a few PS. For t = T/4, 
the relevant eigenvalues of U can be rewritten as 

exp(-i7r/iV2) = ^ (e-*"/^ + exp(-mj;)e*"/^) . 

The corresponding expression for the evolution operator is then 

t/(T/4) = ^ [e"*"/^ -t- e*"/^exp(-mJ^) . 
V2 L J 

For the same initial state |^E'(0)) = |0; J), we apply Eq. (17), and arrive at the result for the 

state at t = T/4 

|^(T/4)) = -^ {e'^^l^ 10; J) + e^-Z^e"- ^ |0 + tt; J)] , (24) 

v2 

which is the linear superposition of two PS's. For times t that are equal to other simple 

rational fractions of the period T it = T j?)^ T /^^..) a superposition of more PS's is found. 
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One may use the fact that the eigenvalues exp(— m/i ^t) of U are periodic in /i with some 
integer period p. Therefore these eigenvalues can be expressed as a finite Fourier series in 
powers of exp(27ri/i/p), which is equivalent to expressing the evolution operator U{t) as a 
finite sum of rotations around the 2;-axis. 

When the number N of particles is odd, so that the values of J and /i are half-integer, 
full revival of the initial state is again found after one period t = T. In fact, since 2/i is an 
odd number, (4/i^ — l)/4 is always an even integer, and it follows that both at time T and 
T/2, the evolution operator is just a phase factor 

U{T) = exp(-i7r/2),t/(T/2) = exp(-i7r/4). 

Hence, apart from a phase factor, full revival is found already at half the time T. In order 
to obtain the evolution operator at the time t = T/4, it is convenient to use the identity for 
half integer values of /i 

exp(-i7r/iV2) = ^e-^"/s (e'"^/^ ^ ^-in,^/2^ ^ 
For the evolution operator this gives the expression 

t/(T/4) = ^e"^"/^ [exp(mJ,/2) + exp(-i7rJ,/2) . 
V2 L J 

For the initial PS |^&(0)) = |0; J), we obtain for the state vector at time T/4 

|^(T/4)) = ^e-'"/s [e'"-'/^ ^^ _ ^/2; J) + e-'"^/^ ^^ ^ ^/2; J)] . 
V2 

Revivals of the state of a BEG have been observed in an optical lattice [13]. 

V. PERIODIC MODULATION OF ENERGY DIFFERENCE 

A simple example of a periodic modulation of the two-well system is to include a time- 
varying energy difference between the two wells. This is realized by substituting in the 
iV-particle Hamiltonian (11) the harmonically varying parameter e {t) = £1 cos cut, while 5 
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and K remain constant. It is convenient to describe the evolution in an interaction picture 
that removes the diagonal terms in the Hamiltonian. We introduce the transformed state 
vector \^'{t)) by the relation 



\^{t))=T{t) M/'(t) 



(25) 



where the state vector I^P (t)) obeys the Schrodinger equation with the Hamiltonian (11), 
and the transformation operator T(t) is defined by 



9i 



T{t) = exp —iO (t) Jz — iKtJ. 



(26) 



with 6{t) = f^dt e (t) = ei (sinut) joj. Notice that the transformed state |^'(t)) has the 
same distribution over the number states \^ as the actual state \^{i)). The transformed 
Schrodinger equation has the standard form 



ih- 



dt 



H'{t) M/' (t) 



(27) 



An explicit form of the transformed Hamiltonian 



H\t) = -h5f\t)J^f{t) 



follows from the general transformation rule [14] 



/(J.)J+ = J+/(J. + 1). 



(28) 



(29) 



Ths relation (29) holds for any analytical function / of the operator J^. After substituting 
(29) into (28), we arrive at the result 



H'{t) = -M 
^ ' 2 



J gi0(t)+«t(2J,+l) ^ Y{,c. 



After a Fourier expansion of the exponentials, we find 



M 



H\t) = -^Y,Jn (^i/^) J+e*n42^^+i)+H + H.c.) . 



(30) 

The form (30) of the operator H'{t) allows a clear physical interpretation. The oscillating 
energy difference e{t) is equivalent to a series of harmonic couplings between the wells with 
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equally spaced driving frequencies nu. The amplitude for each harmonic is proportional 
to the Bessel function of the corresponding order. So the effective coupling between the 
number states |/i) depends strongly on the frequency. 

The Hamiltonian H'{t) contains only non- vanishing elements coupling neighboring num- 
ber states l/i) and |/i + 1). A resonance occurs for the nth harmonic when 

nu + fi;(2/i + 1) = 0, (31) 

which requires that k (2/i + 1) /cj is an integer. 

The strength of this coupling is — i7^J„(£i/eti;)/2, with 

% = (5v/(J-/i)(J + /i+l). (32) 

The effective coupling by the nth harmonic is measured by the parameter 

K = -r—^k—^Jn f ^) , (33) 



'^ ncu + /t(2/i+l) " Vcu/ ' 

which is the ratio of the coupling strength and the detuning from resonance for the transition. 
Whenever |f/^| <^ 1, the coupling is weak. 

When the oscillation frequency u is large compared with the maximal diagonal frequency 
splitting K,{2J + 1), all the time-dependent couplings are weak, and the dominant coupling 
term is the static one with n = 0. The effect of the modulated energy difference is then that 
the coupling term is reduced by the factor Jq^ei/uj). In the high-frequency limit u; ^ ei, 
this factor is one, and we recover the case of a static and symmetric double-well potential 
with El = 0. 

A simple isolated resonance between two number states can occur involving the states 
l/i) with /i = —J or /i = J, since these can be coupled to only one other state. Suppose that 
at t = all atoms are in one of the two wells, so that 

\^it = 0)) = \-J). (34) 
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2 46 8 10 s* 12 24 6 8 

FIG. 1. Time dependence of the populations P^ for the state |/i) = \—J) and |/Li) 
The parameters are taken as 6/ k = 0.25, Ei/k = 14, lo/ n = 3, N = 2J = 16. 



-J + l>. 



If one chooses the frequency u such that the resonance condition (31) holds at certain 
integer no, the corresponding harmonic can be made dominant. Indeed, provided that 
\UJJi\ -C 1 for /i = — J+ 1, for all n, coupling to other states is weak, and we have an effective 
two-level system. This is possible provided that at resonance u is large compared with k, 
which in turn is large compared with the coupling parameter 6. This is demonstrated in 
Fig. 1, where oscillations between the states \—J) and \ — J + 1) are displayed for the initial 
state (34). This means that a single atom out of A^ atoms resonantly oscillates between the 
wells. Upon decreasing the coupling between the wells, the rate of off-resonant coupling 
is decreasing, so one approaches ideal Rabi oscillations between resonant levels. Weaker 
coupling implies a larger oscillation period. The two-level behavior can only occur for 
a system with a nonlinear term J^, since for a linear system the various transitions are 
simultaneously in resonance [6] and [7]. 

In the case that the modulation frequency u is of the same order as k, resonances on 
the different transitions can coincide, and the initial state (34) can spread out over many 
number states. For example, in the simple case that uj = n, the resonance condition (31) 
shows that for each value of //, there is a harmonic n = — (2/i + 1) that is resonant, and the 
population spreads out over all number states. 
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FIG. 2. Time dependence of the fluctuation \Jz) of operator J^ at resonance (a) u; = k and 
out of resonance (b) lo = 6k. Tlie other parameters are taken the same as in the previous figure. 

The difference with the high-frequency case is demonstrated in Fig. 2, where we plot 
the fluctuations AJf of Jz as a function of time, for the initial condition (34), for tu = 4k 
(a) and u = n (b). In the first case, the fluctuations remain limited. In the second case, a 
resonance occurs on each transition, and AJf continues to increase. Even for a very small 
coupling between wells, resonances designed in such a way can lead to enhancement in the 
tunneling rate. This is close to the experimental situation for the double-well trap presented 
in Ref. [4]. Again, this situation is specific for a system with a non-linear term j| in the 
Hamiltonian, since for a linear system various transitions have the same effective coupling. 
Since the coupling is proportional to J„ {ei/u), a resonant transition can be turned off by 
setting the ratio ei/u equal to a zero of the Bessel function. 
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FIG. 3. Time dependence of the expectation value of Jz operator is plotted at resonance to = k. 
The ratio ei/uj is chosen such that (a) ei/uj = 24.26918 (J(i3(ei/u;) = 0) , (b) ei/uj = 23.2759 
(Jii {ei/lo) = 0) (c) Ei/io = 24.93493 {Jj {ei/to) = 0). The total number of particles is chosen 
N = 16. 

This can be used to restrict the evolution to a limited number of states, thereby locating 
a desired number of particles in one of the wells. We demonstrate this idea in Fig. 3. There 
we start from the same initial condition (34) and see that 11 atoms out of 16 are localized 
in the left well if one chooses the ratio Si/uo such that J-j {ei/uo) = (Fig. 3a). Then, taking 
Jii (si/uj) = 0, or Ji5 (ei/cu) = 0, one can localize 13 or 16 particles in one of the wells (Figs. 
3b and 3c). 

VI. GENERALIZATION TO AN OPTICAL LATTICE 

The discussion of the previous section for two wells with an energy difference can be 
generalized to the case of a multimode system, consisting of a chain of potential wells. As a 
model, we take a BEC in a tilted optical lattice [7]. As usual, we neglect the higher bands 
in the lattice, and we consider only a BEC trapped in the lowest energy band, that roughly 
speaking is composed of the ground states in all the wells [15]. If one takes the Wannier 
states |/) with / = ... — 2, —1,0, 1, 2, . . .as the basis of one-particle states, the Hamiltonian 
in second quantization is a direct generalization of Eq. (6) for two wells, and it takes the 
form 

Hbh = — ^ X] i^hi^i + ajai+i) + ^^ ajajaiai + He (t)^ lajai, (35) 

/ I I 

where ai i a] ) are bosonic annihilation (creation) operators in a single Wannier state, 5 and 

K, are the obvious multimode generalizations of two-mode definitions for the nearest neighbor 

coupling and interaction constant (1, 7), 5 is the energy difference in frequency units between 

neighboring Wannier states, which determines the uniform force. This Hamiltonian defines 

the so-called Bose-Hubbard model. 
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The time evolution in a lattice is governed by the time-dependent Schrodinger equation 

(36) 



«i^w>=5,„|*„)>. 



dt 
The uniform force and the interatomic interaction can be eliminated by the substitution 



\'^BH{t))=TBH{t) ^Bi/W 



with 



TBHit) =exp -2^(t)^/: 



aiai-—2_^alal 



am 



I I 

and 6 [t) = f^ dt e (t) is the area of pulse. The Schrodinger equation for the transformed 
state l"^ BH (^)) follows by using the transformation properties of the annihilation operator 

Tg^{t)aiTBH{t) = ai exp ( —il9{t) — int ld\ai — 1 j j = exp ( —il9{t) — iKtajai ) ai, 

which leads to the identity 



Tljf{t)al^^aiTBHit) = aj^idi exp i6{t) + int ia\^{di+i - a\ai + 1 



We obtain the evolution equation 

-dWBHm 



ih- 



with the effective Hamiltonian 



dt 



H'bh I^W it)) 



(37) 



H 



h6 



BH 



-tEK. 



lO; exp 



i6{i) + iKt ( a[_|_]^a;+i — 'di'di + 1 



H.c. 



(38) 



For the case of a periodically modulated uniform force, described by e (t) = Si cosut, this 
Hamiltonian can be put in the form 



H 



BH 



-T^^J'-Vi 






I n 



UJ 



incut + iKt ( ajj^-ydi^i — ajai + 1 



+ H.C.). (39) 



i\^\ = |. . . , A^„i, A^o, Ni, . . .) where two 



This Hamiltonian couples collective number states 

neighboring wells / and / + 1 have exchanged one particle. The coupling between states with 

A*"; = p, Ni^i = q and Ni = p — 1, Ni^i = g + 1 is resonant for a harmonic n when 
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nuj + K{q-p+l) = 0, (40) 

which is analogous to the resonance condition (31). At small tunneling rate we can exclude 
non-resonant coupling terms while assuming that their effective coupling rate is negligible. 
When the uniform force also contains a constant term, so that 6{t) = 6q + 6i cos cut, we have 
to add a term eot to 6{t), and the resonance condition is modified into 

nuj + €o + K{q-p+l) =0. (41) 

When p — q = 1, this condition is independent of k, and a resonant oscillation can occur 
between states with A*"; — A^;+i = ±1. 

Another interesting case is a Mott insulator state, with the same number of particles 
Nq in each well. Such a state has been predicted in [15] and has been recently experimen- 
tally realized in ( [16]), where one (two) atoms have been put in a single lattice site. So, 
|\1' (t = 0)) = |. . . , Nq, Nq, Nq, ...,). Thls statc is directly coupled to the collective Fock state 
which arises if a boson escapes to a neighbouring well, so it has Nq + 1 atoms in one lattice 
site, and Aq — 1 in the neighboring one. Then the resonant condition is nu + £o + ^^ = 0- 

Just as in the case of two wells, resonances coincide when uj is of the same order as 
K. When UJ = K, there is always a harmonic that resonantly couples neighboring wells. 
In the absence of the constant term Eq, the resonance condition takes the universal form 
n + q — p + 1 =0. So, if in the Mott insulator phase the number fluctuations are suppressed 
between wells, we obtain their increase at resonances. 

VII. PERIODIC MODULATION OF COUPLING 

In this Section we consider the effects of a periodic modulation of the coupling coefficient 
6{t) between the wells. As a simple model, we assume that 6 contains a harmonic component, 
so that 

5{t) = 5q + 5i cos cot. (42) 
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In order that the even state \g) = (|1) + |2))/^/2 is the ground state, we keep S{t) positive 
at all times, and we choose 6i to be smaller than 60. Hence we assume that Sq > 61 > 0. 
Moreover, we take the energy of the two wells to be equal, so that e = 0. The Hamiltonian 
in the form of (11) with the coupling coefficient (42) can be easily implemented in practice. 
It describes a BEC in a two-well configuration with a periodically modulated barrier height. 
Precise calculations of the coupling coefficient are given in [5]. 

Since in the Hamiltonian (11) the term proportional to J^ is periodically modulated, 
we expect that the basis of states l/i)^,, which are eigenstates of the operator J^, is the 
natural basis to describe the evolution. Then it is convenient to describe the Hamiltonian 
in terms of the operators J^ and J^, which are defined in Eq. (19). By using the identities 
Jz = —{Jx + ■Jx)/'^ ^^^ -Jx-Jx + -Jx-Jx = 2[J(J + 1) — J^], we rewrite the iV-particle 
Hamiltonian (11) in the form 

H^ = -h6 it) Jx + ^ [jiJ + 1) - J^) + ^ (J^ + J-') . (43) 

This expression demonstrates that a state |/i)^ is coupled only to its next nearest neighbors 
l/i ± 2)^. The coupling strength is measured by the matrix element 

n; = ^x{f^ + 2\ Jf \li)x = ^V(^ + /^ + l)(^ + /^ + 2)(J-/x-l)(J-/x), (44) 

which depends on the interparticle interaction coefficient k, and the particle number N = 2 J. 
In order to get a closer insight to the role of periodic modulation and its resonances, we 
again eliminate the diagonal part of the Hamiltonian, now with respect to the basis of states 
l/x)^. The time-dependent state is expressed 

|vi/(t)) = 5(t)|vl/"(t)), (45) 

with 



S(t) = exp 



iv{t)Jx- -iKt[j{J+l)-l 



(46) 



and Tj it) = Ldt6(t^ is the integrated coupling coefficient. In order to obtain the 
Schrodinger equation for the transformed state |\l'"(t)), we need the transformation prop- 



erty of the off-diagonal operators J^^ and J^ ^. The transformed state |\l/"(t)) obeys the 
Schrodinger equation with the effective Hamiltonian 

In analogy to Eq. (29), we now apply the general transformation rule 

gfyj:' = J^'giZ + 2), (47) 

for an arbitrary analytical function g oi Jx- After making a Fourier expansion we obtain for 
H"{t) the explicit expression 

#"(t) = ^ f] J„ ('^) [j+ 2g-.t(,5o+2.(j.+l)+n.) ^ j^^^ j (48) 



n=—oo 



The form of the Hamiltonian H"{t) resembles the Hamiltonian H'{t), as specified in Eq. 
(30). In the present case, the basis states are the states |/u)^, which are now coupled by the 
square of the corresponding ladder operator J^ ^. As in (30), the coupling term is a series 
of harmonics with equidistant frequencies nuj, with an amplitude proportional to the Bessel 
function of the corresponding order. 

In the high-frequency limit, when the modulation frequency u is large compared with 
the diagonal frequency splittings of the Hamiltonian (43), the effect of the static term pro- 
portional to Jq{26i/uj) in Eq. (48) will be dominant, and the Hamiltonian will be effectively 
constant. Just as in precious sections, the physical reason is that a rapidly modulated field, 
which has a negligible average pulse area, also has a negligible influence. 

On the other hand, Eq. (48) immediately shows that the coupling between the state |/i)^ 
and l/i + 2)^ of the nth harmonic is resonant when 

nuj + 2fi;(/x + I) + 6o = 0. (49) 

The other coupling terms are negligible when the coupling strength is small compared with 
the oscillation frequency, which leads to the weak-coupling criterion 



^^ . f26^ 



'•Jn 



nuj + 2k {ji + 1) + 5q \ uj 
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< 1. (50) 




FIG. 4. Time dependence of the populations Pj^i of the state \ii)^ with /x = J and ji = J — 2. 
The other parameters are taken as (5i/k = 14, 5q/ k = 16, uo/ k = 20, N = 14. 

So, if the initial state |\I/(0)) = \J)^ is resonantly coupled to |J — 2)^, while further 
couplings of this alter state are negligible, we have an effective two-level system. This 
is demonstrated on the figure 4, where besides the resonant oscillations, one obtains non- 
resonant escape of population to the rest of manifold. Upon decreasing the coupling between 
the state |/u)^, only the resonant states are involved and they exhibit clear Rabi oscillations. 
This shows how resonances lead to an escape of population from the initial state to the 
other states in the manifold of states |/i)^. The similarity with the response to the periodic 
modulation in the form of a periodically modulated energy difference between the wells. 
Recall that the state |J)^ is the state in which all particles are in the even state \g) = 
(|1) + |2))/-\/2, which is the one-particle ground state. In the state \J — 2)^, two particle 
have been transferred to the odd excited state. 

VIII. CONCLUSIONS 

A BEC trapped in a two-well potential can be expected to be very sensitive to the fre- 
quency of any applied periodic perturbation. We test this idea by periodically modulating 
the asymmetry or the barrier height of such a configuration. Compared with the analogous 
situation of a single atom trapped in a light field with a periodic modulation, the many- 
particle nature of the BEC gives rise to some new effects. For both types of modulation, 
two-state resonances may be observed, where a single atom out of the BEC oscillates be- 
tween the wells. It is also possible to enter a regime of parameters where more than two 
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states are resonantly coupled, with more than one particle oscillating between the wells. 
Using such resonances, one can manipulate the average number of particles in the wells by 
varying the relevant parameters, such as the magnitude and the modulation frequency of 
the energy difference. This effect can be considered also as a means to resonantly enhance 
the tunnelling rate between wells. We generalize the basic ideas developed for two-wells to 
a multiwell system, such as a BEC in an optical lattice. Whereas the periodic modulation 
of the energy difference is related to coupling between number states in the two wells, the 
periodic modulation of the height of the barrier leads to coupling between number states in 
superposition states of the two wells with specific values of the relative phase. 
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